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ABSTRACT 

We have observed two blank fields of approximately 30 by 23 arcminutes using the 
William Herschel Telescope. The fields have been studied as part of the Canadian 
Network for Observational Cosmology Field Galaxy Redshift Survey (CNOC2), and 
spectroscopic redshifts are available for 1125 galaxies in the two fields. We measured 
the lensing signal caused by large scale structure, and found that the result is consistent 
with current, more accurate measurements. 

We study the galaxy-galaxy lensing signal of three overlapping samples of lenses 
(one with and two without redshift information), and detect a significant signal in all 
cases. The estimates for the velocity dispersion of an = 0) = 5.6 x 10® /i“^Lb 0 
galaxy agree well for the various samples. The best fit singular isothermal sphere model 
to the ensemble averaged tangential distortion around the galaxies with redshifts yields 
a velocity dispersion of cr* = 130lF km/s, or a circular velocity of V* = 184+25 km/s 
for an Lg galaxy, in good agreement with other studies. 

We use a maximum likelihood analysis, where a parameterized mass model is 
compared to the data, to study the extent of galaxy dark matter halos. Making use 
of all available data, we find cr* = 111 ± 12 km/s (68.3% confidence, marginalised 
over the truncation parameter s) for a truncated isothermal sphere model in which 
all galaxies have the same mass-to-light ratio. The value of the truncation parameter 
s is not constrained that well, and we find s* = 260+/3"' h~^ kpc (68.3% confidence, 
marginalised over cr*), with a 99.7% confidence lower limit of 80 k~^ kpc. Interestingly, 
our results provide a 95% confidence upper limit of 556 k~^ kpc. The galaxy-galaxy 
lensing analysis allows us to estimate the average mass-to-light ratio of the field, 
which can be used to estimate flm- The current result, however, depends strongly on 
the assumed scaling relation for s. 

Subject headings: cosmology: observations — dark matter — gravitational lensing 


1 INTRODUCTION 

The (small) differential deflection of light rays by interven¬ 
ing structures allows us to study the projected mass distribu¬ 
tion of the deflectors, without having to rely on assumptions 
about the state or nature of the deflecting matter. The first 
attempt to detect this effect, called weak gravitational lens¬ 
ing, was made by Tyson et al. (1984), who tried to measure 
the signal induced by an ensemble of galaxies. This area of 
astronomy blossomed with the successful measurements of 
the signal induced by rich clusters of galaxies at interme¬ 
diate redshifts (e.g., Tyson, Wenk, & Valdes 1990; Bonnet, 
Mellier, & Fort 1994; Fahlman et al. 1994; Squires et al. 

* Based on observations made with the William Herschel Tele¬ 
scope operated on the island of La Palma by the Isaac Newton 
Group in the Spanish Observatorio del Roque de los Muchachos 
of the Institute de Astrofisica de Canarias. 


1996; Luppino & Kaiser 1997; Hoekstra et al. 1998; for an 
extensive review see Mellier 1999). 

These studies of rich clusters were an important first 
step in demonstrating the feasibility of weak lensing anal¬ 
yses, but nowadays more and more studies concentrate on 
blank fields. For example, galaxy groups have masses inter¬ 
mediate between clusters of galaxies and galaxies. Hoekstra 
et al. (2001) measured the ensemble averaged weak lensing 
signal from a sample of 50 groups identified by Carlberg et 
al. (2001) in the Canadian Network for Observational Cos¬ 
mology Field Galaxy Redshift Survey (CNOC2). 

Other applications of wide field lensing are the mea¬ 
surement of the lensing signal caused by large scale struc¬ 
ture (Bacon et al. 2000,2002; Hoekstra et al. 2002a, 2002b; 
Kaiser, Wilson, & Luppino 2000; Refregier et al. 2002; van 
Waerbeke et al. 2000, 2001, 2002; Wittman et al. 2000), 
and the study of galaxy biasing (Hoekstra, Yee & Gladders 
2001b; Hoekstra et al. 2002c). Another important applica- 
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tion is the study of the dark matter halos of field galax¬ 
ies (e.g., Brainerd, Blandford, & Smail 1996; Griffiths et al. 
1996; Dell’Antonio & Tyson 1996; Hudson et al. 1998; Fis¬ 
cher et al. 2000; Wilson, Kaiser, & Luppino 2001; McKay et 
al. 2001; Smith et al. 2001). 

Rotation curves of spiral galaxies have provided impor¬ 
tant evidence for the existence of dark matter halos (e.g., 
van Albada & Sancisi 1986). Also strong lensing studies of 
multiple imaged systems require massive halos to explain the 
oberved image separations. However, both methods provide 
mainly constraints on the halo properties at relatively small 
radii. The weak lensing signal can be measured out to large 
projected distances, and in principle it can be a powerful 
probe of the potential at large radii, constraining the extent 
of the dark matter halos (e.g., Brainerd et al. 1996, Hudson 
et al. 1998; Fischer et al. 2000). Only satellite galaxies (e.g., 
Zaritsky & White 1994) provide another way to probe the 
outskirts of isolated galaxy halos. 

The lensing signal induced by an individual galaxy is 
too low to be detected, and one has to study the ensemble 
averaged signal around a large number of lenses. Redshifts 
for the individual galaxies are useful, because they allow a 
proper scaling of the lensing signal around the galaxies, and 
they are necessary for studies of the evolution of the mass-to- 
light ratio of field galaxies from lensing. Hudson et al. (1998) 
were the first to make use of (photometric) redshifts in their 
galaxy-galaxy lensing analysis of the northern Hubble Deep 
Field. Unfortunately, the small area covered by the HDF 
limited the accuracy of their results. 

The analysis of commissioning data of the Sloan Digital 
Sky Survey (SDSS) by Fischer et al. (2000) was a major 
step forward. Fischer et al. (2000) detected a very significant 
lensing signal, demonstrating the importance of the survey 
for the study of galaxy halos. More recently, McKay et al. 
(2001) used the available redshift information from the SDSS 
to study the galaxy-galaxy lensing signal as a function of 
galaxy properties. 

We obtained deep R-band imaging data for two fields 
that have been studied as part of the CNOC2 Field Galaxy 
Redshift Survey (e.g., Yee et al. 2000). Earlier results on 
groups of galaxies, based on these data, were presented by 
Hoekstra et al. (2001a). In this paper, we use the data to 
study the galaxy-galaxy lensing signal of three overlapping 
samples of galaxies (one with, and two without redshift in¬ 
formation) . 

The structure of the paper is as follows. In Section 2 we 
present the observations and data reduction. In this section 
we also describe in detail the object analysis and the correc¬ 
tions for the various observational distortions. In Section 3 
we discuss the redshift distribution of the sources we use in 
this study. We investigate the lensing by large scale structure 
in Section 4. The analysis of the galaxy-galaxy lensing signal 
is presented in Section 5. In Section 6 we present our esti¬ 
mates of the field mass-to-light ratio for different halo mod¬ 
els. Throughout the paper we take Hq = lOOh km/s/Mpc, 
Dm = 0.2, and Da = 0, although the results do not depend 
critically on the adopted cosmology. 


field 

RA 

DEC 

date 

1447 

14'*47™04.3'’ 

09°12'39" 

May 19 - 21 1998 

2148 

21'* 48”* 34'* 

-05°56'00" 

Aug 30 - Sep 2 1997 


Table 1. Positions of the centres of the fields used for the weak 
lensing analysis. The last column gives the dates of the observa¬ 
tions. 


2 DATA 

The Canadian Network for Observational Cosmology Field 
Galaxy Redshift Survey (CNOC2) targeted four widely sep¬ 
arated patches on the sky to study the field population 
of galaxies in the universe. Redshifts of ~ 6200 galaxies 
with a nominal limit of Rc = 21.5 were measured, result¬ 
ing in a large sample of galaxies at intermediate redshifts 
{z = 0.12 — 0.55). A detailed description of the survey, 
and the corresponding data reduction is given in Yee et al. 
(2000). In this paper we study the dark matter properties of 
the galaxies targeted by CNOC2, and to this end we make 
extensive use of the redshifts and multi-colour photometry 
obtained by the CNOC2 survey. 

We observed the central parts of the two CNOC2 
patches 1447+09 and 2148-05 (Yee et al. 2002, in prepa¬ 
ration) using the 4.2m William Herschel Telescope (WHT) 
at La Palma. The images were taken using the prime focus 
camera, equipped with a thinned 2048 x 4096 pixels EEVIO 
chip, with a pixel scale of 0'.'237 pixeU^. The resulting field 
of view of the camera is approximately 811 by 16l2. 

The patches observed in the CNOC2 survey are much 
larger than the field of view of the WHT prime focus camera, 
and we observed a mosaic of 6 pointings. Table 1 lists the 
central positions of the observed fields as well as the dates of 
the observations. The typical integration time per pointing 
is one hour in R (see Table 2). 

2.1 Data reduction 

The images were flatfielded, using a master flatfield con¬ 
structed from the science exposures. The images were cal¬ 
ibrated using observations of standard stars from Landolt 
(1992). 

The data for each pointing typically consists of 3 expo¬ 
sures of 1200s, which were taken with small offsets. Table 2 
lists the total integration time, seeing, and number of de¬ 
tected objects for both observed fields. Note that because 
of adopted weighting scheme for the weak lensing analysis, 
the “effective” number of galaxies is approximately 30% of 
the total number. The data for the 1447 field have a me¬ 
dian seeing of ~ 0"7, whereas the data for the 2148 field are 
somewhat worse, with a seeing of 0”85. 

The exposures had to be remapped before the images 
were combined into the final image, because of focal plane 
field distortions. We selected and measured the positions of 
stars in each of the exposures and used these as input for 
the IRAF tasks geomap and geotran. To obtain the final 
images that were used for the object analysis, the remapped 
images were simply averaged to ensure that neither cosmic 
ray rejection or medianing changed the shape of the PSF or 
the galaxies in a non-linear way. 

Although a remapping of the images was necessary, the 
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field 

(1) 

^exp 

[s] 

(2) 

seeing 

[”] 

(3) 

^ objects 

(4) 

n 

[arcmin”^] 

(5) 

mim 

(90%) 

(6) 

1447-1 

6000 

0.75 

5738 

46 

24.3 

1447-2 

3600 

0.70 

5545 

44 

24.5 

1447-3 

3600 

0.65 

6419 

51 

24.8 

1447-4 

3600 

0.70 

5570 

45 

24.6 

1447-5 

4800 

0.70 

6087 

49 

24.5 

1447-6 

3600 

0.65 

6826 

55 

24.8 

2148-1 

3600 

1.0 

4712 

38 

24.2 

2148-2 

3600 

0.80 

4976 

40 

24.3 

2148-3 

3600 

0.75 

5540 

44 

24.2 

2148-4 

4800 

0.85 

4804 

38 

24.2 

2148-5 

3600 

0.90 

4486 

36 

24.1 

2148-6 

3600 

0.85 

5110 

41 

24.2 


Table 2. Summary of the deep WHT imaging. Column (1): iden¬ 
tification of the pointing; (2) total integration time; (3) median 
seeing; (4) number of galaxies; (5) number density of galaxies; (6) 
90% completeness limit in R. 

camera induced distortion is small. The WHT prime focus 
observer’s manual lists the coefficents to estimate the tele¬ 
scope distortion, which is purely radial. We also calculated 
the distortion from a comparison of positions of bright stars 
that coincide with stars from the USNO catalog. The rel¬ 
atively large uncertainty (~ (X'S) in the astrometry of the 
USNO catalog limits the accuracy of this approach. The re¬ 
sults, however, agree with the distortion derived from the 
parameters listed in the WHT observer’s manual, and we 
use the latter to calculate the camera distortion used in our 
analysis. We find that the induced shear is small: at most 
0.66% in the corners of the field. The corrected weak lensing 
distortion field is obtained by subtracting the camera dis¬ 
tortion from the observed distortions (which are corrected 
for PSF anisotropy), as described in Hoekstra et al. (1998). 

2.2 Object detection and analysis 

Our analysis techique is based on that developed by Kaiser, 
Squires, & Broadhurst (1995) and Luppino & Kaiser (1997), 
with a number of modifications which are described in Hoek¬ 
stra et al. (1998). We analyse each pointing separately as 
the seeing and PSF anisotropy vary between each exposure. 
After the catalogs have been corrected for the various ob¬ 
servational effects, they are combined into a master catalog 
which covers the complete field that is observed. 

The first step in the analysis is to detect the faint galaxy 
images, for which we used the hierarchical peak finding al¬ 
gorithm from Kaiser et al. (1995). The peak finder also pro¬ 
vides an estimate for the Gaussian scale length Vg of the 
object. We select objects which were detected with a signif¬ 
icance V > 5a over the local sky in the combined images. 
This galaxies are used for the weak lensing analysis. We also 
run the peak finder on the images of the single exposures. 
We do a coincidence test on these catalogs and classify ex¬ 
tremely small, but very significant objects as cosmic rays, 
which are removed from the catalog that is used for the ob¬ 
ject analysis. 

Some faint cosmic rays may hit galaxies, and conse¬ 


quently might not be recognized as cosmic rays. Based on 
the number of cosmic ray hits, and the area covered by galax¬ 
ies we find that less than 0.2% of the galaxies might be af¬ 
fected. Also, cosmic rays only introduce additional noise in 
the shape measurement, but do not bias the result. Con¬ 
sequently we conclude that remaining cosmic rays have a 
negligble effect on our results. 

For the weak lensing analyis we select objects which are 
detected in at least two of the shorter exposures. The objects 
which are detected in only one of the shorter exposures are 
small faint objects, which are not useful as their shape pa¬ 
rameters are noisy. The resulting catalogs are inspected vi¬ 
sually in order to remove spurious detections, such as spikes 
from saturated stars, HIT regions in resolved galaxies, etc. 

For all detected objects we measure the apparent mag¬ 
nitude in an aperture with a radius of 3 times the Gaus¬ 
sian scale length of the object, the half light radius, and 
the shape parameters (polarization and polarizabilities). We 
also estimate the error on the polarization following Hoek¬ 
stra, Franx, & Kuijken (2000). 


2.3 PSF correction 

To measure the small, lensing induced, distortions in the im¬ 
ages of the faint galaxies it is important to correct the shapes 
for observational effects, such as PSF anisotropy and seeing: 
PSF anistropy can mimic a lensing signal, and the correc¬ 
tion for seeing is required to relate the measured shapes to 
the real lensing signal. 

To do so, we follow the procedure outlined in Hoekstra 
et al. (1998) (which is based on Kaiser et al. 1995; Luppino 
& Kaiser 1997). We select a sample of moderately bright 
stars from our observations. These are used to character¬ 
ize the PSF anisotropy and seeing. Figure 1 shows a typical 
result for one of the pointings. We fit a second order polyno¬ 
mial to the observed shape parameters of the stars, and this 
model is used to correct for the PSF anisotropy. Figure lb 
also shows the residual polarization of the stars after the 
correction. The residuals are small indicating that we can 
reliably correct for the PSF anisotropy. 

The next step is to correct the shapes for the circular¬ 
ization by the PSF. The stars that were used to study the 
PSF anisotropy are also used to compute the ‘pre-seeing’ 
shear polarizability (Luppino & Kaiser 1997; Hoekstra 
et al. 1998). The measurement of P^ is very noisy, and we 
combine the estimates of many galaxies to reduce the noise. 
We also found that the size of the PSF, and thus the cor¬ 
rection, depends on the position on the chip. The variation, 
however, is small: about 10% maximum. To account for this, 
we bin the raw polarizabilities not only in bins of Vg, but also 
as a function of position. For a given rg we fit a second order 
polynomial to the median P'^’s as a function of position, and 
use the results to compute the for each galaxy. Figure 2a 
shows P^ as a function of rg. 

Objects with small values for P^ require large correc¬ 
tions, thus increasing the noise. The weighting scheme sug¬ 
gested in Hoekstra et al. (2000) gives already less weight 
to these objects, but in addition we exclude objects with 
P^ < 0.1 
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Figure 1. Shape of the PSF in one of the pointings of the 1447 field; (a) the polarizations of the stars before correction; (b) residual 
polarization after correcting for the PSF anisotropy using a second order polynomial fit to the data, (c) PSF anisotropy as a function 
of position on the chip. The left panel shows the PSF anisotropy before correction. To show the higher order components the average 
polarization has been subtracted. The direction of the sticks indicates the major axis of the PSF, whereas the length corresponds to the 
size of the anisotropy. The right panel shows the residuals in the stellar polarizations after correction. 


3 REDSHIFTS OF LENSES AND SOURCES 

For a given mass distribution, the amplitude of the weak 
lensing signal is proportional to the dimensionless mass sur¬ 
face density 


•^crit 

where the critical surface density is defined as 

““ 47rGAAs 47rGA/3’ 

Here Ds, D\, and D\s correspond to the angular diameter 
distances between the observer and the source, observer and 
the lens, and the lens and the source. The lensing signal 
depends on both the redshifts of the lenses and the sources, 
and the dependence on the source redshift is characterized 
by P, which is defined as 

P = max[0, As/A]. (3) 

If the redshift of the lens approaches that of the source, 
the amplitude of the lensing signal, which is proportional to 
P, decreases. 

Most of the lensing signal comes from galaxies that are 
generally too faint to be included in redshift surveys. To re¬ 


late the lensing signal to a physical mass it is necessary to 
know the redshift distribution of the faint background galax¬ 
ies. To this end, we use photometric redshift distributions 
from the Hubble Deep Fields North and South (Fernandez- 
Soto, Lanzetta, & Yahil 1999). Hoekstra et al. (2000) used 
these results for their analysis of the distant cluster MS 1054- 
03 {z = 0.83), and concluded that these redshift distribu¬ 
tions provide a good approximation of the true distribution. 
We use the colours of the galaxies in the HDFs to derive 
their R band magnitude. For each lens-source pair we com¬ 
pute the corresponding value of p. 

The lenses considered in our analysis are at much lower 
redshifts than MS 1054-03, and consequently the uncertainty 
in the value of P is small. Based on fleld-to-field variation in 
the redshift distribution and the uncertainty due to the finite 
number of galaxies in the Hubble Deep Fields, we estimate 
that the uncertainty in P is 2%. 


4 LENSING BY LARGE SGALE STRUGTURE 

One of the selection criteria for the CNOC2 fields (which are 
described in Yee et al. 2000) is that no known nearby rich 
cluster should be in the field. However, the observed fields 
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5000 

§ of galaxies 


Figure 2. (a) as a function of rg for all galaxies. Due to seeing 
variations the curves for the various fields blend; (b) number of 
galaxies with a given (c) number of galaxies with a given 
value of P'l'. For the weak lensing analysis we restrict the sample 
to galaxies with P'^ > 0.1. 


might contain distant massive clusters. In principle, such 
clusters can be found in a weak lensing analysis, provided 
they are massive enough (e.g., Wittman et al. 2001), even 
if they are “dark” (see Erben et al. 2000). We note that 
projection effects can actually introduce spurious detections 
(e.g., Hoekstra 2001; White, van Waerbeke & Mackey 2002). 

To investigate possible structures in the CNOC2 fields, 
we reconstructed maps of the dimensionless surface density, 
using the original Kaiser & Squires (1993) algorithm. The 
resulting mass maps are consistent with noise maps. We find 
a few 3 c 7 peaks, but no obvious counterparts are seen in the 
number counts of bright galaxies. The mass maps have been 
smoothed with a Gaussian with a FWHM of 1 arcminute, 
resulting in approximately 500 independent points in each 
map. Therefore one expects about three 3 ct peaks. 

It is useful to estimate the mass detection limit for our 
data. We assume that the cluster mass profile is well de¬ 
scribed by a singular isothermal sphere (SIS): 


where ve is the Einstein radius. The observed scatter in the 
shapes of the sources results in a typical uncertainty in the 
Einstein radius r_B of ~ 2". We fit SIS models out to 5 ar- 
cminutes, but note that the uncertainty in the determination 
of the Einstein radius does not depend significantly on this 
particular range. Hence, we should be able to detect a clus¬ 
ter with te > 6” (at the 3 (t level). The velocity dispersion 
of the SIS model is given in km/s by 


a = 186.3 



km/s. 


(5) 


when r_B is given in arcseconds. Given the redshift distribu- 


field 

(91 > 

{ 92 } 

^9 

1447 

-0.00334 

0.00144 

0.00226 

2148 

0.00323 

-0.00140 

0.00254 


Table 3. Average distortion in the observed fields. The errors in 
the measurements are given in the last column. The values are 
small, which suggest that the correction for the PSF anisotropy 
has worked well. 



Figure 3. The observed variance of ( 7 ^) as a function of the 
radius of the aperture in which the shear is averaged. Note 
that the points are correlated. The error bars do not include 
the additional uncertainty caused by cosmic variance. For ref¬ 
erence we have displayed the expected variances for three differ¬ 
ent cosmologies: SCDM {flm = 1.0; Oa = 0; F = 0.5; crs = 0.5), 
OCDM (n^ = 0.3; Qa = 0; F = 0.21; cts = 0.85), and ACDM 
(n^ = 0.3; flA = 0.7; F = 0.21; as = 0.9) 


tion of our sources a cluster with a velocity dispersion of 630 
km/s at z = 0.2 would be detectable, whereas a cluster at 
2 : = 0.5 would have to have a velocity dispersion larger than 
900 km/s. Such rich clusters would have been detected in 
the CNOC2 redshift survey. A detectable cluster at 2 > 0.5 
would have to be even more massive. Such massive systems 
are very rare, and are unlikely to be in the observed fields. 
Based on the mass reconstructions we find no massive clus¬ 
ters in the CNOC2 fields. 

Although we do not detect significant structures in the 
GNOC2 field, any large scale structure along the line of sight 
will introduce an excess alignment in the shapes of the faint 
galaxies (compared to a random field). The measurement of 
this signal provides a direct measurement of the statistical 
properties of the large scale mass distribution (e.g., Bland- 
ford et al 1991; Kaiser 1992; Bernardeau, van Waerbeke, & 
Mellier 1997; Schneider et al. 1998). 

A widely used method is to look for the excess variance 
in an aperture: the top-hat variance. This measurement as 
a function of scale can be compared to predictions from cos- 
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mological models. The measurement is difficult, and resid¬ 
ual systematics (such as an imperfect correction for PSF 
anisotropy) increase the observed variance. 

Since the first detections (Bacon et al. 2000; Kaiser, 
Wilson, & Luppino 2000; Maoli et al. 2001; van Waerbeke 
et al. 2000; Wittman et al. 2000) tremendous progress has 
been made. The most recent results are based on large data 
sets, and yield consistent results, despite the small statistical 
errors (e.g.. Bacon et al. 2002; Hoekstra et al. 2002a, 2002b; 
Refregier et al. 2002; van Waerbeke et al. 2001, 2002). The 
area covered by the two CNOC2 fields is smaller than the 
area covered by these studies, but it is still interesting to 
examine our results and compare them to the predictions: 
although the analysis will not provide tight constraints on 
the cosmology, it is a useful test whether residual systematics 
are present in the data. 

For the analysis, we use galaxies with apparent mag¬ 
nitudes 21 < i? < 26 as our sample of sources. The mea¬ 
surements of the average distortion of the two fields are pre¬ 
sented in Table 3. To measure the excess top-hat variance 
as a function of aperture radius, we tile the observed fields 
with apertures of a given scale, compute the variance, and 
subtract the contribution from the intrinsic shapes of the 
galaxies (e.g.. Bacon et al. 2000; Kaiser et al. 2000; Hoek¬ 
stra et al. 2002a; van Waerbeke et al. 2000, 2001). 

The results are presented in Figure 3. The error only 
includes the statistical uncertainty caused by the intrinsic 
shapes of the sources, and does not include the contribution 
from cosmic variance. Also note that the points are corre¬ 
lated. The measured variances are small, and consistent with 
current, more accurate measurements. 

The top-hat variance is very sensitive to residual PSF 
anisotropy, and can be used to test the accuracy of the cor¬ 
rection. For instance, if we do not correct the shapes of the 
galaxies for PSF anisotropy, we measure an excess variance 
( 7 ^) = (8.4 ± 0.6) X 10“® at a scale of 1 arcminute, 60 times 
larger than the signal presented in Figure 3. If we subtract 
90% (under-correction) or 110% (over-correction) of the PSF 
anisotropy, the observed variance increases by a factor ~ 2 . 
Because any residual PSF anisotropy will increase the sig¬ 
nal, this test indicates that the adopted correction yields the 
minimum (or close to) variance. 


5 GALAXY-GALAXY LENSING 

We use two different methods to study the lensing signal 
caused by the field galaxies. First we measure the ensem¬ 
ble averaged tangential distortion around the lens galaxies 
(galaxy-mass correlation function), which provides a robust 
estimate of the lensing signal. We also perform a maximum 
likelihood analysis, which is described in Section 5.3. 

It is important to note that the measurement of the 
galaxy-galaxy lensing signal is much less sensitive to a small 
residual systematics: we measure the lensing signal that is 
perpendicular to the line connecting many lens-source pairs. 
These connecting lines are randomly oriented with respect to 
the PSF anisotropy, and hence suppress any residual system¬ 
atics. Given the results obtained for the top-hat smoothed 
variance (see Section 4), we expect that the effect of residual 
systematics is negligible. 

In Table 4 we present the three subsamples of lenses 


that we study. The ‘faint’ and the ‘bright’ sample are magni¬ 
tude limited, and we use statistical redshift distributions for 
the lenses to infer the average halo properties. The ‘CNOC2’ 
sample has the same magnitude limits as the ‘bright’ sam¬ 
ple, but consists only of galaxies with spectroscopic red- 
shifts. It includes approximately half of the galaxies of the 
‘bright’ sample. The ‘faint’ sample is comparable to the sam¬ 
ple of lenses studied by Brainerd et al. (1996) who used 
20 < riens < 23, and 23 < rsource < 24. Brainerd et al. 
(1996) did not use fainter sources, because they did not cor¬ 
rect for observational distortions. 


5.1 Lenses selected irrespective of redshift 
information 

The ensemble averaged tangential distortion as a function 
of radius from the lens is a well established way to present 
the galaxy-galaxy lensing signal (e.g., Brainerd et al. 1996; 
Hudson et al. 1998; Fischer et al. 2000; McKay et al. 2001). 
The results for the ‘faint’ and the ‘bright’ sample are shown 
in Figure 4. For both samples a significant lensing signal is 
detected. If the measured signal is caused by gravitational 
lensing, no signal should be observed when the sources are 
rotated by 7r/4. The results of this test are shown in figure 4b 
and d, and no signal is seen indeed. 

If the lenses are distributed randomly on the sky, the 
tangential distortion profile can be related directly to the 
ensemble averaged mass profile of the lenses. In reality the 
lenses cluster, and the observed signal is the convolution 
of the mass profile and the galaxy correlation function: 
the tangential shear profile measures the galaxy-mass cross¬ 
correlation function. As a result mass estimates based on the 
tangential shear profile are biased high, because at large dis¬ 
tances from the lens, one measures the mass of the lens and 
associated galaxies. However, the lens dominates the signal 
on small scales. 

We fit a SIS model to the tangential distortion pro¬ 
files presented in figure 4. We do not know the redshifts of 
the individual lenses in these samples, and consequently we 
can only determine the ensemble averaged Einstein radius 
(r_B). For the ‘faint’ sample we find a best fit Einstein radius 
{te) = (/.'118 ± (y.'025, and for the ‘bright’ sample we obtain 
(rs) =0'.'176±0'.'045. 

We examined the effect of an imperfect correction for 
PSE anisotropy on the determination of the Einstein ra¬ 
dius. Even in the extreme case that no correction is applied 
the derived value is changed by only 5%. In the case of a 
90% or 110% correction of the PSF anisotropy, the signal 
is changed by ~ 1%. Based on the results obtained for the 
top-hat smoothed variance (see Section 4), we conclude that 
the effect of an imperfect correction of the PSF anisotropy 
is much smaller than 1%, and hence is negligible. 

We use the effective /? for these samples (see column 
6 of Table 4) to derive a mass weighted average velocity 
dispersion = 1081^2 km/s for the ‘faint’ sample, 

and = 122lj® km/s for the ‘bright’ sample. The 

corresponding circular velocity can be calculated using 14 = 
ypia. 

The derived values of depend on the selection of 

the sample of lens galaxies, and one cannot compare these 
results to findings of other studies, given the differences in 
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Figure 4. (a) The ensemble averaged tangential distortion for the ‘faint’ sample. The solid line corresponds to the best SIS model fit, 
which has te = 0^118 di 0^^025. Panel c shows the signal around the ‘bright’ sample, and we find ve = 0^176 di 07045. Panels b and d 
show the average signal when the sources are rotated by 7r/4. No signal should be present if the signal in panels a and c is caused by 
lensing. 
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[”] 

[”] 

[km / s] 

[km/s] 

[km/s] 

faint 

20 < R < 23 

23 < R < 26 

8715 

0.46 

0.35 ± 0.01 

0.118 ± 0.026 

0.159 ± 0.035 

126+]3 

1781]] 

1241]] 

bright 

17.5 < R< 21.5 

22 < R < 26 

2125 

0.34 

0.41 ± 0.01 

0.176 ± 0.046 

0.161 ± 0.042 

1171]® 

1651]] 

1151]] 

CNOC2 

17.5 < R< 21.5 

22 < R < 26 

1125 

0.36 

0.40 ± 0.01 

0.196 ± 0.047 

0.196 ± 0.048 

130l]® 

1841]] 

1291]] 


Table 4. Properties and results for the different samples of lens galaxies. (2) the range in apparent magnitude of the lens galaxies; (3) 
range in apparent magnitude for the sources; (4) number of lens galaxies; (5) median redshift of the lenses; (6) average value of /3 based 
on the redshift distributions of the lenses and the sources. (7) best fit Einstein radius; (8) estimate for the Einstein radius of an Lg 
galaxy; (9) best estimate for the velocity dispersion of an Lg galaxy under the assumption that the luminosity evolves oc (1 + z)\ (10) 
corresponding circular velocity; (11) the velocity dispersion of an Lg galaxy for no luminosity evolution. 


sample selection. Instead we estimate the velocity dispersion 
(or circular velocity) of an Lg galaxy. We assume a scaling 
relation between the velocity dispersion and the luminosity 
of the galaxy of the form 

G oc ^ or Vc oc . ( 6 ) 

We also assume that the luminosity of a lens of given 
mass evolves with redshift ex (1 + 2 :) (e.g., Lin et al. 1999). 
With these assumptions the average value of the Einstein 
radius (in radians) is given by 


{te) 


4 ^ /o.fW\ 

= 0) \^V l + «/’ 


(7) 


where is the observed intrinsic luminosity of the galaxy. 
We also introduce {r*g) 


^e) 



( 8 ) 


which is the Einstein radius of an Lg galaxy at the average 
redshift of the sample of lenses. 

The redshift distributions of the sources and the ‘faint’ 
lens sample are derived from the photometric redshift dis¬ 
tributions of the HDF North and South ((Fernandez-Soto, 
Lanzetta, & Yahil 1999). For the ‘bright’ sample we use 
the redshift distribution from the CNOC2 survey with the 
proper weighting to take into account the incompleteness of 
the survey. To determine the restframe B luminosities we 
use template spectra for a range in spectral types and com¬ 
pute the corresponding passband corrections as a function 
of redshift and galaxy colour (this procedure is similar to 
the one described in van Dokkum & Franx 1996). 

Lin et al. (1999) have studied the field galaxy luminos¬ 
ity function at intermediate redshift from the CNOC2 sur¬ 
vey, and derived Mg{z = 0.3) = —19.18 -I- 51ogh, which 
corresponds to a luminosity of 7.3 x 10® h^^Le©. With 
our assumed luminosity evolution, this results in Lg (2 = 
0) = 5.6 X 10® /i“®Lb 0 . Madgwick et al. (2002) derived 
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Mg = —19.55 + 5 log ft from the the 2dF redshift survey, 
which probes lower redshifts. Hence, the CNOC2 value (Lin 
et al. 1999) is rather low, in particular with our choice of lu¬ 
minosity evolution. Throughout the paper we assume that 
the luminosity evolves oc {1 z), but for reference we also 

list the results for the no-evolution scenario, in which case 
Lb(2 = 0) = 7.3 X 10® ft-^Le© . 

Column 8 in Table 4 lists the derived Einstein radii 
of an Lb galaxy for the ‘bright’ and the ‘faint’ sample. 
For the ‘bright’ sample the observed {ve) is actually very 
close to r*E- For the ‘faint’ sample we find that the observed 
{ve) = 0.74(rJ;). Columns 9 and 10 list the velocity disper¬ 
sion and circular velocity of an L^ galaxy under the assump¬ 
tion that the luminosity evolves oc {1 + z). Column 11 shows 
the resulting velocity dispersion without luminosity evolu¬ 
tion. Coincidentally, the values for a* for the evolving and 
non-evolving case are very similar. Furthermore, the results 
for the faint and bright sample agree well with one another. 


5.2 Lenses with redshifts from CNOC2 

Redshifts for the lens galaxies are useful because they al¬ 
low a proper scaling of the signals around the lens galax¬ 
ies. Intrinsically faint (and therefore low mass) galaxies, or 
galaxies with redshifts comparable to the source galaxies are 
given lower weights. We scale the observed distortion of each 
galaxy, as well as its error, such that it corresponds to that of 
an Lb galaxy at the median redshift of the ‘CNOC2’ sample 

...I. ^'5.6xlO"(l + ^)V'Vo.‘l^ .1. 

=(— r, - ) 

where we assumed that the luminosity scales with the fourth 
power of the velocity dispersion. The value of /3 is calculated 
for each galaxy separately, based on its redshift, and the 
redshift distribution of the sources. The sample of lenses 
consists of 1125 galaxies with spectroscopic redshifts, and 
17.5 < Riens < 21.5. We note that, because of the selection 
of these galaxies, the sampling in apparent magnitude is not 
uniform (the completeness for the fainter galaxies is lower). 
Figure 5a shows the resulting average tangential distortion 
as a function of radius around an L^ galaxy at a = 0.36. A 
significant galaxy-galaxy lensing signal is detected. 

The best fit SIS model to the scaled tangential distor¬ 
tion yields an Einstein radius te = 0Vl96±0V047. From this 
result we derive an average velocity dispersion of cr = 130!ti7 
km/s or circular velocity of Vc — 184^25 km/s for an L^ 
galaxy. If we assume no luminosity evolution, we obtain 
cr* = 129^^7 km/s (with Lb (2 = 0) = 7.3 x 10® ft^^Ls© ). 

The number of lenses in the ‘CNOC2’ sample is half 
that of the ‘bright’ sample, but the signal-to-noise ratio with 
which the lensing signal is detected is similar. Thus redshifts 
for the lenses are useful as they can be used to improve 
the signal-to-noise ratio of the measurement. A comparison 
of the results for the different samples (which are listed in 
Table 4) shows a good agreement. 

We note that we have made assumptions about the lu¬ 
minosity evolution and the scaling relations. The redshift 
information for the lenses is crucial when one would like to 
constrain the scaling relations or study the evolution of the 
mass-to-light ratio of these field galaxies. 


radius [h * Mpc] 



radius [”] 


Figure 5. (a) The ensemble averaged tangential distortion 
around galaxies with 17.5 < R < 21.5 and measured spectro¬ 
scopic redshifts, scaled to an Lg galaxy at a redshift z = 0.36. 
The solid line corresponds to the best SIS model fit, which has 
Te = Ol'lOO ± CK'047. (b) The average signal when the sources 
are rotated by n/A. No signal should be present if the signal in 
panel a is caused by lensing. 


5.3 Sizes of galaxy halos 


As mentioned above, the clustering of lenses will bias the 
mass estimates inferred from the observed tangential distor¬ 
tion profile. In this section we use a parameterised model 
for the mass distribution of individual galaxies and compare 
the predicted distortions to the observed distortion field. We 
use a truncated halo model, in order to constrain both the 
mass and the extent of the dark matter halos. 

This approach naturally takes into account the cluster¬ 
ing of the lenses in the sample. Furthermore, we make use of 
both components of the distortion. It is important to note 
that the contribution of associated faint galaxies is not ac¬ 
counted for. The effect of a large-scale cosmological shear on 
the parameter estimation is negligble because one measures 
the lensing signal for many randomly oriented lens-source 
pairs. The effect of clusters/groups on the results is exam¬ 
ined in Section 5.4. 

A useful model to describe a truncated halo is (Schnei¬ 
der & Rix 1997) 


E(r) 


2Gr 



r 

\/r^ + 


( 10 ) 


where s is a measure of the truncation radius. The total 
mass of this model is finite, and half of the mass is contained 
within r = |s. The total mass is given by 


Mtot = = 7.3 X 10^^ 

Gr 


100 km/s J \ 1 Mpc 


( 11 ) 


The scatter in the polarizations of the galaxies is ap¬ 
proximately constant with apparent magnitude, and it can 
be well approximated by a Gaussian distribution. In that 
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case the log-likelihood follows a distribution with the 
number of degrees of freedom equal to the number of free 
model parameters (e.g., Hudson et al. 1998), and the de¬ 
termination of confidence intervals is straightforward. The 
log-likelihood is given by the sum over the two components 
of the polarization a of all the source galaxies 


log r = - 


E 


5 * (^* ) ® * ) P? 


( 12 ) 


where gij are the model distortions, is the shear polar¬ 
izability, Cij are the image polarizations for the jth galaxy, 
and (Je^ is the uncertainty in the shape measurement of the 
jth galaxy. 

In our maximum likelihood analysis we ignore the con¬ 
tribution from lenses outside the held of view (e.g., Hudson 
et al. 1998). For small helds of view this tends to lower the 
resulting a* and s* somewhat. The area covered by our ob¬ 
servations is much larger than the HDF North studied by 
Hudson et al. (1998), and we hnd that the effect on our 
estimates of ct* and s* is negligible. 

To scale the signals of the various lenses we use again 
a oc which is based on both dynamical and observa¬ 

tional considerations. The situation is different for the trun¬ 
cation parameter s, because there are no observational con¬ 
straints. Here we will explore several cases. 

If all halos have the same value for s the total mass- 
to-light ratio scales as (M/L)tot oc (where we as¬ 

sume that L oc a^). Thus the mass-to-light ratios of more 
luminous galaxies are lower. Another option is to take 
(M/L)tot = constant for all galaxies. This choice is equiv¬ 
alent to taking s oc (e-g-, Brainerd et al. 1996; Hudson 
et al. 1998). The last relation we examine is s oc cr'^, which 
gives (M/L)tot oc L1/2. 

Figure 6 shows the results for the ‘CNOC2’ sample, for 
which the individual redshifts of the lenses are known. It 
shows the likelihood contours for the parameters ct* and s* 
jointly. The result for a constant s for all galaxies is pre¬ 
sented in figure 6a, and the result for s oc is shown in hg- 
ure 6b. We omit the likelihood plot s oc a^, but we list the 
best fit parameters (68.3% confidence limits) for all three 
models in Table 5a. The results change slightly when we 
assume no luminosity evolution, and the values of the pa¬ 
rameters are listed in Table 5d. The conhdence intervals on 
(7* (s*) listed in Tables5a and d are obtained by marginalis¬ 
ing over St (n*). Note that these confidence intervals cannot 
be inferred directly from Figure 6 (e.g.. Press et al. 1992, 
Numerical Recipes in C, Figure 15.6.3). 

We also examined the influence of an imperfect correc¬ 
tion for PSF anisotropy. The change in the derived velocity 
dispersion is approximately 1% if we under or overcorrect 
by 10%. The value of s* is changed by ~ 6% in these case. 
As we have demonstrated above, the actual correction for 
PSF anisotropy is much better than 10%, and hence we hnd 
it to be negligible. 

In all cases the velocity dispersion is well constrained, 
and in good agreement with the results from the ensemble 
averaged tangential distortion. The value of the truncation 
parameter has signihcant freedom, and the value of s* de¬ 
pends strongly on the assumed scaling relation. We hnd that 
the value of s* decreases when s scales with a high power of 
the velocity dispersion. The minimum for fhe three mod¬ 


els are comparable, and no scaling relation for s is preferred 
over the other ones. Table 5a also lists the best estimate and 
the 68.3% conhdence limits for the total mass and mass-to- 
light ratio of an Lggalaxy. The total galaxy mass is not well 
constrained, mainly because of the large uncertainty in the 
value of the truncation parameter s*. 

The mass model uses only half of the galaxies with 
17.5 < R < 21.5, and therefore it ignores the contribu¬ 
tions of the other galaxies (they are given zero mass). If 
the remaining galaxies were distributed randomly, their con¬ 
tribution to the lensing signal is that of noise. We know, 
however, that galaxies are clustered. Consequently the mass 
model will associate the mass of the galaxies without red- 
shifts with the galaxies in the ‘CNOC2’ sample that have 
measured redshifts. 

To examine this in more detail we study the ‘bright’ 
sample, which is selected irrespective of redshift informa¬ 
tion. Even this selection is not ideal, neighbouring galaxies 
that are fainter than the applied magnitude limits are still 
excluded. If a substantial fraction of the mass is in these 
galaxies both the values of cr* and s* are overestimated. 

Figure 7 shows the resulting likelihood contours. To ob¬ 
tain figure 7 we have assumed that all halos have the same 
velocity dispersion, and the same value for s in arcseconds. 
We hnd = 117ljg km/s, and (s) = 881*^27 arcseconds 

(68.3% conhdence). At the average redshift of the lenses, this 
corresponds to (s) = 284jtgy®h“^kpc. The estimate for (s) 
is a weighted average over the value of s for all the galaxies, 
which are all at diherent redshifts. The interpretation of the 
result is not straightforward for the various options for the 
scaling relation of s. 

To infer the best estimates for cr*, and s* of the bright 
sample one has to perform a maximum likelihood analysis 
in which the redshift of each individual galaxy is a free pa¬ 
rameter, which has to be chosen such that it maximizes the 
likelihood. This approach is computationally not feasible, 
and we use another approach to obtain estimates for cr* and 
s*. 

We use the same method that was used for the galaxies 
with redshifts, but instead we create mock redshift catalogs. 
The redshift distribution of the galaxies in the ’bright’ sam¬ 
ple has been measured by the CNOC2 survey. The redshifts 
of the galaxies in the mock catalogs are drawn randomly 
for the CNOC2 survey based on the apparent R magnitude 
of the galaxies in the ’bright’ sample. We take the incom¬ 
pleteness of the survey into account when the redshifts are 
drawn from the survey. This mock redshift catalog allows 
us to determine the best estimates for cr* and s*, based on 
a maximum likelihood analysis of the observed distortion 
held. We repeat this procedure 25 times, and use the av¬ 
erage x^ surface. The best estimates for cr* and s* for the 
three different scaling relations for s are listed in Table 5b. 

We verihed that this procedure yields an unbiased es¬ 
timate, by simulating the data and applying the procedure 
described above. We used the derived values of cr* and s* 
(see Table 5b) as input values for the simulations. These 
simulations also allow us to estimate the scatter in the mea¬ 
surements. Compared to the situation where the redshift of 
each lens is known, the lack of redshifts increases the uncer¬ 
tainty in St by ~ 50%. The uncertainty in ct* is increased 
by 30%. For instance, adopting s oc cr^, the uncertainty in- 
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Figure 6. Likelihood contours for the velocity dispersion cr* of an L*g(z = 0) = 5.6 X 10® h ‘^LgQ galaxy, and its truncation radius s*. 
The contours indicate 68.3%, 95.4%, and 99.7% confidence limits on two parameters jointly. The cross indicates the best parameters. 
We use O' oc and two different scaling relations for s. The dashed lines give models with a total mass of 1, 5,10 X IO^^/i^^Mq. (a) 

The truncation radius s is the same for all galaxies, (b) The results when s oc cr^, which corresponds to the case where all galaxy have 
the same total mass-to-light ratio. 


creases by ~ 45h~^ kpc and the uncertainty in u* increases 
by ~ 3 km/s. 

Comparison of the results preseuted in Table 5b with 
the results obtained for the ‘CNOC2’ sample (Table 5a) 
shows that the estimates for ct, and s* are smaller, but 
consistent with one another. For the case that s oc 
the ’bright’ sample yields cr* = 111 ± 15 km/s and s* = 
220/lioo^~^ kpc, compared to (t, = 126tig km/s and s* = 
313/lgj®/i“^ kpc for the ’CNOC2’ sample. The agreement 
with the results from the direct averaging method is good. 

The analysis of the ‘bright’ sample has ignored any 
available redshift information. To make optimal use of the 
available information, we redo the analysis as follows. For 
the galaxies with redshifts we use the observed values, but 
for the remaining galaxies we draw redshifts and luminosities 
from the CNOC2 survey randomly (as was done previously 
for all galaxies in the ‘bright’ sample). As before, we create 
mock catalogs which are analysed. This procedure improves 
the accuracy in the estimates of a, and s*. The results are 
presented in Table 5c. Simulations show that the lack of 
redshifts for all galaxies with 19.5 < 7? < 21 from CNOC2 
increases the uncertainty in a* by 18% and the uncertainty 
in St by 10%, compared to the situation where full redshift 
information is available. 


5.4 Effect of galaxy groups 

It is well known that galaxies cluster, and most of the galax¬ 
ies reside in groups of galaxies. If the matter in galaxy groups 
is associated with the halos of the group members (i.e. these 
halos are indistinguishable from the halos of isolated galax¬ 
ies) the analysis presented above gives a fair estimate of the 
sizes of galaxy halos. However, if a significant fraction of 


the dark matter in galaxy groups is distributed in a com¬ 
mon group halo, the interpretation of the results becomes 
difficult. 

Fischer et al. (2000) argued that galaxy groups compli¬ 
cate the interpretation of the galaxy-galaxy lensing signal. 
They measured the lensing signal out to large projected dis¬ 
tances, and find that the distortion does not vanish at large 
radii, which is caused by large scale clustering (the angular 
correlation function declines relatively slowly). 

Compared to Fischer et al. (2000) we have the advan¬ 
tage that a number of galaxy groups have been identified in 
the CNOC2 fields (Carlberg et al. 2001). The weak lensing 
signal of these groups was studied by Hoekstra et al. (2001a). 
They derived a mass weighted group velocity dispersion of 
km/s. 

We will use the groups to examine their effect on the 
galaxy-galaxy lensing results for the galaxies with redshifts 
from the CNOC2 survey. In addition to the galaxies, we in¬ 
clude the lensing signal from group halos to our mass model. 
The groups are placed at their observed positions, and the 
groups are modeled as singular isothermal spheres. To study 
the effect of the groups on the galaxy-galaxy lensing results 
we increase the velocity dispersion of the ‘groups’ from 0 
km/s (the result without groups) up to 300 km/s. For each 
choice of the group velocity dispersion we perform a max¬ 
imum likelihood analysis and determine the best estimates 
for CT* and s*. 

For (Tgroup < 150 km/s the best estimates for cr* and 
St vary by a few percent. For larger group velocity disper¬ 
sions the value for cr* and s* decrease slowly with increas¬ 
ing CTgroup. For a group velocity dispersion of 300 km/s, the 
value of a* has dropped by ~ 10%, and the value of s* has 
decreased by ~ 20%. We note that the minimum has 
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Table 5. (a) Results from the maximum likelihood analysis using the galaxies with redshifts from the CNOC2 survey, (b) The results for 
the analysis of the ‘bright’ sample, (c) Results when we use the observed redshifts if available, and the redshift distribution if the redshift 
has not been measured, (d)-(f) list the results for no luminosity evolution, for which case Lg (2 = 0) = 7.3 x 10^ /i~^Lb 0 • Only the 
results presented in (a) and (d) are real maximum likelihood parameters. See the text for more details. The errors correspond to 68.3% 
confidence. (1) the scaling relation used for the truncation parameter s; (2) estimate for cr*; (3) estimate for s*; (4) 99.7% confidence 
lower limit on s*; (5) 95% confidence upper limit on s*; ( 6 ) estimate for the total mass; (7) total mass-to-light ratio of an Lg galaxy. 


increased significantly for cTgroup = 300 km/s, compared to 

O-group = 0 km/s. 

Ideally one would like to use different halo models for 
the group members and the ‘isolated’ galaxies and study 
the difference in the best parameters for these two types of 
galaxies. With the current data we cannot perform such an 
analysis, because the number of group members is too low. 
However, our results are based on approximately one quarter 
of the full CNOC2 data set. An analysis of the full survey 
will improve the signal-to-noise ratio of the measurements 
by a factor ~ 2, and will allow a better study of the effect 
of galaxy groups. Also numerical simulations can be useful 
to examine if, and how one can separate the contribution of 
the galaxy halos and the smooth group halos. 


5.5 Comparison with other lensing studies 

Other studies selected different samples of lenses, used dif¬ 
ferent scaling relations, or made other assumptions about 
the luminisoty evolution. Fortunately, sufficient information 
is available to allow for a useful, if crude, comparison. We 
will compare the findings of some other studies to the results 
from our maximum likelihood analysis. 

The typical lens galaxy in the sample studied by Hud¬ 
son et al. (1998) is at 2 ; ~ 0.6, and has Mb = -18.5-1-5 log h 
(for qo = 0.5). Hudson et al. (1998) derived a circular veloc¬ 
ity of 210 ± 40 km/s. This corresponds to a galaxy with a 
luminosity oi Lb{z = 0) ~ 3.2 x 10® h~^LBQ for our choice 
of cosmology and luminosity evolution. Our weak lensing 
analysis suggests a circular velocity of 14 = 119 ± 12 km/s 
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Figure 7. Likelihood contours for the velocity dispersion cr, of an 
Lg{z = 0) = 5.6 X 10® galaxy, and the average value of 

the truncation radius {s) a determined from the ‘bright’ sample. 
We have also indicated the physical scale when s is the same 
for all galaxies. The contours indicate 68.3%, 95.4%, and 99.7% 
conhdence limits on two parameters jointly. The cross indicates 
the best parameters. The dashed lines give models with a total 
mass of 1, 5,10 X lO^®/i“^M 0 . 

for such a galaxy. If we assume no luminosity evolution (as 
Hndson et al. (1998) did), we obtain a circular velocity of 
Vc = 130 ± 14. The two results are inconsistent at the ~ 2a 
level. It is not clear, however, how to interpret the differ¬ 
ence in circular velocity, as our weak leasing analysis probes 
a larger physical scale than Hudson et al. (1998), and the 
mix of galaxy types might also be different. 

Fischer et al. (2000) measured the galaxy-galaxy lensing 
signal with much higher accuracy than Hudson et al. (1998) 
from SDSS data. Fischer et al. (2000) find for their sam¬ 
ple of lenses an average mass weighted velocity dispersion 
of = 145 — 195 km/s (95% confidence). McKay et 

al. (2001) noted that the SDSS redshift survey showed that 
the lens redshifts used by Fischer et al. (2000) were overesti¬ 
mated by ~ 35%. The correct range in line-of-sight velocity 
dispersion of the Fischer et al. (2000) analysis should be 
« 105 - 145 km/s. 

Most of the signal comes from galaxies with luminosities 
around Lgi = 8.7 x lO®/i“^Lg/0 (P. Fischer, private com¬ 
munication). If we use an average B — V = 0.55 and the 
transformations from Fukugita et al. (1996), we find that 
the adopted B band luminosity of our Lg galaxy corre¬ 
sponds to a g' band luminosity of Lg/ = 6 x lO®/i“^Lg/0. 
Our lensing analysis implies a velocity dispersion of ct = 
122jjJ^2 km/s (68.3% confidence) for a galaxy with Lgi = 
8.7 X lO®h“®Lg/0, which agrees well with the Fischer et al. 
(2000) result. 

McKay et al. (2001) used a sample of galaxies with spec¬ 
troscopic redshifts from the SDSS. Their sample of lenses is 
comparable in size to the one used by Fischer et al. (2000). 
Using the complete sample, McKay et al. (2001) find a ve¬ 


locity dispersion a = 100 — 130 km/s (95% confidence) 
with a best fit value of cr = 113 km/s for a galaxy with 
Lg/0 ~ 9 X lO®h.“^Lg/0. For such a galaxy, our lensing anal¬ 
ysis yields a = 122/1^2 km/s (68.3% confidence), in excellent 
agreement with McKay et al. (2001). 

Fischer et al. (2000) and McKay et al. (2001) also at¬ 
tempted to constrain the sizes of galaxy halos. Both stud¬ 
ies indicate that the the galaxy halos are large. Unfortu¬ 
nately it is not clear what the correct value for the Fischer 
et al. (2000) result is. McKay et al. (2001) used the same 
approach as Fischer et al. (2000) and find a lower limit of 
^mm _ 230 /i“^ kpc (95% confidence). McKay et al. (2001) 
do not use any scaling for s (they assume all halos are the 
same) and their results should be compared to our results 
for constant s for all galaxies. We obtain an estimate of 
s = s* = 337 ^/qq/i“^ kpc, which is in good agreement with 
McKay et al. (2001). 

Interestingly, McKay et al. (2001) find that the mass 
within an aperture of radius 260h“^ kpc scales proportional 
to the luminosity of the lens. This appears to be in contra¬ 
diction with the naive expectation from the Tully-Fischer or 
Faber-Jackson relation, which suggests M oc \/L. However, 
the latter is no longer true for truncated halos. As discussed 
in Section 5.3, the relation between the mass and luminosity 
on large scales depends on the adopted scaling relation for 
s (at large radii M oc a^s). For instance, s (x a'^ results in 
M oc L at radii (much) larger than s. 

Hence, the results presented by McKay et al. (2001) 
suggest that the truncation is smaller than 260h~^ kpc. This 
is likely to be the case for the faint lenses, but the brighter 
lenses are expected to have larger halos. However, luminous 
galaxies are clustered more strongly (e.g., Norberg et al. 
2002) and the contribution from neighboring galaxies could 
bias the lensing signal somewhat high, although the effect is 
expected to be small (McKay et al. 2001). 

Recently, Yang et al. (2002) used numerical simulations 
to study the results presented by McKay et al. (2001). The 
galaxies in the simulations obey the Tully-Fischer or Faber- 
Jackson relation. Yang et al. (2002) find that the aperture 
mass within 260h,“^ kpc is proportional to the luminosity 
of the lens. The results from Yang et al. (2002), however, 
suggest that the linear relation found by McKay et al. (2001) 
is a coincidence, because they find that the halo mass scales 
oc ® (which corresponds to s oc ct^). The mass with an 
aperture of radius 260h~^ kpc results in a linear relation, 
whereas other apertures would yield different results. 

Finally we note that the results presented here are in 
good agreement with preliminary results from an analysis of 
R band imaging data from the Red-Sequence Cluster Survey 
(RCS). Hoekstra, Yee, & Gladders (2001d) use lenses with 
19.5 < R < 21. Their results imply cr* = 102 ± 4 km/s and 
s* = 198ljg/i~^ kpc for an Lg(z = 0) = 5.6 x 10® h“^LB0 
galaxy (adopting s oc cr^), which is in good agreement with 
our results. 

5.6 Comparison with Tully-Fisher relation 

Comparison of the luminosity and the rotation velocity of 
gas in spiral galaxies has shown that there exists a tight re¬ 
lation between these two observables: the Tully-Fisher (TF) 
relation. Note, however, that weak lensing probes the mass 
on a much larger scale than the rotation curve. For example. 
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if the rotation curve is slowly declining, the lensing result 
will be lower (compared to the flat rotation curve). Further¬ 
more, the sample of lenses used hereis different from the 
sample of spiral galaxies used for the TF relation. Ideally 
one would like to select a sample of spiral galaxies for this 
comparison, but this is outside the scope of this paper. Nev¬ 
ertheless, it is interesting to compare our lensing result to a 
local determination of the TF relation. 

One of the most recent determinations of the TF rela¬ 
tion was presented by Verheijen (2001), who studied a large 
sample of galaxies in the Ursa Major cluster of galaxies. The 
sample has the advantage that the galaxies are at a distance 
of 18.6 Mpc (Tully & Pierce 2000). This allowed Verheijen 
(2001) to study the scatter in the TF relation. 

Verheijen (2001) has measured Vfiat, the rotation veloc¬ 
ity of the ‘flat’ part of the rotation curve, which is repre¬ 
sentative for the mass of the galaxy. As is customary in TF 
studies, Verheijen (2001) used extinction corrected luminosi¬ 
ties. The luminosities used in our lensing analyses, however, 
have not been corrected for inclination and extinction. To 
allow a direct comparison with the weak lensing results, we 
use the uncorrected magnitudes from Tully et al. (1996). 
The resulting TF relation is presented in Figure 8 (open 
dots with error bars). 

The hatched region indicates the Icr interval around the 
TF relation, based on the scatter in the observations. In the 
comparison we use the circular velocity corresponding to the 
value of O', (assuming s oc o^) listed in Table 5c. We assume 
that I4iat = The result is indicated by the large black 
dot. The comparison depends on the value of the Hubble 
parameter, because the inferred luminosities of the lenses 
in our analysis depend on h, whereas the absolute magni¬ 
tudes of the galaxies from Verheijen (2001) are based on the 
adopted distance to the Ursa Major cluster (and therefore 
do not depend on the Hubble parameter). For the compar¬ 
ison presented in Figure 8 we adopted a value of Ho = 70 
km/s/Mpc. Although one cannot make a direct comparison 
between the weak lensing results and the TF relation (for 
reasons mentioned above) it is comforting that the results 
are rather similar. 


6 MASS-TO-LIGHT RATIO AND Um 

A useful method to estimate the matter density of the 
universe was proposed by Oort (1958); Qm is the prod¬ 
uct of the average mass-to-light ratio of the universe, and 
its luminosity density. Carlberg et al. (1997) measured the 
mass-to-light ratios of rich clusters of galaxies, and inferred 
Um = 0.19 ± 0.06. The galaxy properties of rich clusters 
are different from those of the field, and a large correction is 
needed to relate the cluster mass-to-light ratio to that of the 
field. Smaller corrections are required when galaxy groups 
are used. Hoekstra et al. (2001a) derived = 0.19 ± 0.1 
from their weak lensing analysis of 50 groups. Bahcall & 
Comerford (2002) found Q,m = 0.17 ± 0.05 from a combina¬ 
tion of clusters and groups. 

A measurement of the average mass-to-light ratio of 
field galaxies has the advantage that no additional correc¬ 
tions for the difference in stellar populations are required. 
An estimate for the average mass-to-light ratio of the field, 
i.e., the universe as a whole, can be obtained from the re¬ 



Figure 8. Plot of rotation velocity of the ‘fiat’ part of the ro¬ 
tation curve Vflat versus the absolute B-band magnitude (uncor¬ 
rected for inclination and extinction). The open dots with error 
bars indicate the measurements of Verheijen (2001), using the 
(uncorrected) magnitudes from Tully et al. (1996). The hatched 
region indicates the la interval around the derived TF relation. 
In the comparison we use the circular velocity corresponding to 
the value of cr* (assuming s oc a'^) listed in Table 5c. The result is 
indicated by the large black dot. The absolute magnitudes of the 
galaxies from Verheijen (2001) are based on the adopted distance 
to the Ursa Major cluster, and therefore do not depend on the 
value of the Hubble parameter. The inferred luminosities of the 
galaxies studied in the weak lensing analysis do depend on Hq, 
and for the comparison we take a value of Hq = 70 km/s/Mpc. 
We have also indicated how the weak lensing result would change 
with h. The weak lensing result is in good agreement with the 
observed TF relation. We note that weak lensing probes the mass 
on a much larger scale than the rotation curve, and one therefore 
might expect only a qualitative agreement. Furthermore, the sam¬ 
ple of lenses is rather different from the sample of spiral galaxies 
used for the TF relation. 


suits of our maximum likelihood analysis The inferred total 
mass-to-light ratio, however, depends on the assumed scal¬ 
ing relation for the truncation parameter s (as is apparent 
from column 7 in Table 5). 

Lin et al. (1999) have determined the luminosity func¬ 
tion of field galaxies from the CNOC2 survey, and we use 
their results to estimate the average mass-to-light ratio of 
the field. The results are listed in Table 6 (column 4). When 
s (X a^, the average mass-to-light ratio equals that of an Lg 
galaxy. In the case that all halos have the same value for s, 
we already found that the resulting mass-to-light ratio of an 
L* galaxy is high, but in addition, fainter galaxies have even 
higher mass-to-light ratios. This results in a very high field 
mass-to-light ratio. On the other hand, if s oc faint galax¬ 
ies have low mass-to-light ratios, and the field mass-to-light 
ratio is somewhat lower than that of an galaxy. 

To estimate the luminosity density of the universe we 
use the results from Lin et al. (1999), which are also based 
on the CNOC2 survey. We convolve the redshift distri- 
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Table 6 . Estimates for the average field mass-to-light ratio for different scaling relations of the truncation parameter s. The errors 
correspond to 68.3% confidence. (1) the scaling relation used for the truncation parameter s; (2) the dependence of the total mass-to- 
light ratio on the luminosity; (3) total mass-to-light ratio of an Lg galaxy; (4) average field mass-to-light ratio at z = O', (5) average field 
mass-to-light ratio at the average redshift of the sample of galaxies {z = 0.34); ( 6 ) estimate for the matter density of the universe Qm- 


bution of the galaxies with the redshift dependent lumi¬ 
nosity density from Lin et al. (1999), which yields j = 
(3.0±0.6) X lO^hLsoMpc”®. We use this result, and our es¬ 
timates for the average field mass-to-light ratio at 2 = 0.34 
to derive the corresponding values of Qrn, which are listed 
in column 6 of Table 6. 

With the current data, the statistical error in the value 
of flm is large. In addition, the uncertainty in the scaling 
relation of s limits the determination of flm through this 
technique. We also assumed that all matter in the universe 
is associated with the galaxy dark matter halos. The leasing 
signal is not changed by adding a sheet of constant surface 
density (Gorenstein, Shapiro, & Falco 1988). Thus a uni¬ 
formly distributed form of dark matter cannot be detected 
in our analysis. Consequently, our estimate for flm should be 
interpreted as a lower limit. Another approach is to use in¬ 
dependent measurements of flm (e.g., from CMB or cosmic 
shear studies) to constrain the scaling relation for s. 


7 CONCLUSIONS 

We observed two blank fields of approximately 30 by 23 ar- 
cminutes using the William Herschel Telescope. The fields 
were studied as part of the Canadian Network for Observa¬ 
tional Cosmology Field Galaxy Redshift Survey (CNOC2) 
(e.g., Yee et al. 2000; Carlberg et al. 2001), and spectro¬ 
scopic redshifts are available for 1125 galaxies in the two 
fields. Earlier results on groups of galaxies, based on these 
data, were presented by Hoekstra et al. (2000a). 

We have measured the lensing signal caused by large 
scale structure. We observed signal is low, and consistent 
with more accurate measurements from much larger surveys 
(e.g.. Bacon et al. 2002; Hoekstra et al. 2002a,b; Refregier 
et al. 2002; van Waerbeke et al. 2002). 

We examined the effect of an imperfect correction for 
PSF anisotropy, and find the correction used here has 
worked well: it gives the smallest cosmic shear signal. In 
addition we find that the uncertainty in the correction has 
a negligible effect on the galaxy-galaxy lensing results. 

We have studied the ensemble averaged tangential dis¬ 
tortion (galaxy-mass correlation function) around three sub¬ 
samples of lenses. In all three cases a clear lensing signal is 
detected. We relate the lensing signal to an estimate of the 
velocity dispersion of an Lg galaxy, and find that the re¬ 
sults for the three samples agree well with each other. For 
the sample of galaxies with redshifts from the CNOC2 sur¬ 
vey we derive a, = 13041® km/s (or K* = ISdIas km/s). 


Note, however, that this measurement is slightly biased , 
because of the clustering of the lenses. 

To study the properties of the dark matter halos sur¬ 
rounding the lens galaxies, we used a maximum likelihood 
analysis. This technique allowed us to derive constraints on 
both the velocity dispersion and the extent of the dark mat¬ 
ter halos (it naturally accounts for the clustering of the 
lenses). The value of the truncation parameter s depends 
strongly on the assumed scaling relation, but the velocity 
dispersion is well constrained. Galaxy groups can compli¬ 
cate the interpretation of the truncation parameter s. We 
examined how galaxy groups affect our results, and we find 
that their effect is relatively small. Our data are not suf¬ 
ficient for a detailed study, but the full CNOC2 data set 
can be useful to this end, as well as studies of numerical 
simulations. 

Under the assumption that all galaxies have the same 
total mass-to-light ratio, we find a value of cr, = 111 ± 
12 km/s (68% confidence, marginalised over the truncation 
parameter s) for the velocity dispersion of an Lg galaxy. 
For the truncation parameter we obtain s* = 260473“^ 
kpc (68% confidence, marginalised over a*), with a 99.7% 
confidence lower limit of 80 h~^ kpc, and a 95% confidence 
upper limit of 556 h~^ kpc. For this model we find that 
the average field mass-to-light ratio at 2 = 0 is 3934i24 
Mq/Lb 0 , comparable to what is found for rich clusters of 
galaxies (Carlberg et al. 1997) and galaxy groups (Hoekstra 
et al. 2001a). 

The field mass-to-light ratio provides an estimate of Q,m, 
the matter density of the universe. The results from the 
weak lensing analysis provide a lower limit on Qrn, because 
lensing cannot detect dark matter that is distributed uni¬ 
formly through the universe. Unfortunately, our estimate of 
the field mass-to-light ratio, and consequently the derived 
value of Qrn, depends strongly on the assumed scaling of s. 

The results presented here demonstrate that weak lens¬ 
ing is a powerful tool to study the dark matter halos of 
field galaxies, as it can probe the mass distribution out to 
large projected distances. In this paper we have used ap¬ 
proximately one quarter of the CNOC2 survey. An analysis 
of the full survey will improve the accuracy of the results 
by a factor ~ 2. Furthermore, many large imaging surveys 
are currently underway. First results from the Red-Sequence 
Cluster Survey (RCS) indicate that these surveys can pro¬ 
vide accurate constraints on the dark matter halos of field 
galaxies (e.g., Hoekstra et al. 2002d). 
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